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Abstract 

We extract the x dependence of the valence nonsinglet u — d distribution 
function in the nucleon from the lowest few moments calculated on the lattice, 
using an extrapolation formula which ensures the correct behavior in the 
chiral and heavy quark limits. We discuss the implications for the quark 
mass dependence of meson masses lying on the J PC = 1 Regge trajectory. 
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I. INTRODUCTION 



Parton distribution functions (PDFs) contain a wealth of information on the nonper- 
turbative structure of the nucleon. Quark and gluon distributions probed in deep inelastic 
scattering and other high energy processes have provided valuable insights into the workings 
of QCD in the low energy domain. The observation of an asymmetry between d and u quarks 
in the proton sea to take just one example, has served to highlight the important role 
that the pion cloud of the nucleon and the chiral symmetry of QCD [[J plays in hadronic 
structure, even at high energies. 

More generally, studies of PDFs can help with the task of identifying the appropriate 
effective degrees of freedom of QCD at low energies. Through the application of the operator 
product expansion (OPE) to QCD, high energy processes such as deep inelastic scattering 
can be factorized into short and long distance contributions, allowing one to calculate the 
former in perturbation theory, while isolating all of the nonperturbative physics in the 
latter. Over the past two decades considerable experience has been accumulated with various 
nonperturbative, low energy models of the nucleon which have been used to study PDFs 
Initial studies focused on the valence quark distributions as a means of constraining 
valence quark model parameters, although recently more ambitious efforts have attempted 
to describe sea quark and gluon distributions from low energy models. 

Although the model studies have been helpful in exploring the relationship between 
high energy processes and low energy phenomenology, ultimately one would like a more 
exact connection of PDFs with QCD. A mathematically more rigorous approach is provided 
through lattice QCD. Indeed, the determination of the moments of the PDFs is one of the 
benchmark calculations of hadron structure in lattice QCD. Modern computational advances 
have allowed large scale simulations to be undertaken which are progressively improving the 
errors associated with finite lattice spacings and finite volume effects. However, until recently 
H large differences between lattice results and experiment have remained. 

Because PDFs are light cone correlation functions, it is not possible to calculate them 
directly on the lattice in Euclidean space. Instead one calculates moments of PDFs, defined 
(for Bjorken x) as: 

(x n ) g = jf 1 dx x n (q(x) + (-l) n+1 q(x)) , (1) 

which are related through the OPE to matrix elements of local twist-two operators. A num- 
ber of calculations of PDF moments have been performed over the last decade, most notably 
by the QCDSF group || in the quenched approximation. More recently, the MIT group 
jTJ has confirmed the earlier quenched results, and in addition made the first unquenched 
simulations. These results indicate that at the relatively large quark masses at which the 
calculations were made, the unquenched results are indistinguishable from the quenched 
within the current errors. 

Despite the impressive progress of lattice calculations of moments of PDFs, there has 
been a long standing problem in reconciling the lattice data with experiment, which has 
posed a serious threat to the credibility of current lattice calculations. Namely, for the un- 
polarized moments all of the calculations to date, which have been made at quark masses 
of between 30 and 190 MeV, have yielded results which are typically 50% larger than the 



2 



experimental values when linearly extrapolated to the physical quark masses. This discrep- 
ancy was recently resolved with the observation [||] that a linear extrapolation in quark mass 
omits crucial physics associated with the long range structure of the nucleon in the form of 
the pion cloud. In particular, the inclusion of the nonlinear, nonanalytic dependence on the 
quark mass in the extrapolation of the moments from the region of large quark masses to 
the physical, light quark masses allows one to fit both the lattice data and the experimental 
values with good accuracy ||. It therefore restores confidence in the extraction of all hadron 
observables from current simulations. 

In this paper we go one step further and ask what can the lattice data tell us about the x 
dependence of the PDFs themselves? In particular, we shall extract the valence nonsinglet 
distribution u v — d v by extrapolating the lattice data to the physical region using a form 
that ensures the correct behavior in the chiral and heavy quark limits. In Section II we 
examine the initial question of the extent to which the x dependence can be determined 
from just a few low moments of a PDF. After establishing the accuracy to which the PDF 
can be reconstructed, we then extract the valence u v — d v distribution from the calculated 
lattice moments in Section III. Using the chiral extrapolation formula, we compute the x 
dependence of the u v — d v distribution as a function of quark mass. In Section IV we discuss 
the implications of the results for masses of mesons lying on the J PC = 1 Regge trajectory. 
Finally, in Section V we make some concluding remarks. 



II. RECONSTRUCTION OF DISTRIBUTIONS FROM MOMENTS 

The reconstruction of the complete x dependence of a PDF in principle requires infinitely 
many moments, so that a priori the task of building up a PDF from only a few, low moments 
may seem rather ambitious. This is especially true for PDFs at large x, information on which 
is contained in higher moments. On the lattice, because of problems with operator mixing 
for operators with spin > 5, all lattice calculations have so far been restricted to n < 3. It 
is crucial, therefore, if one is to ever make use of lattice data to learn about PDFs, to see 
which features of the PDFs can be reconstructed from just the lowest 3 or 4 moments. 

Formally, the x dependence of the PDF can be obtained by performing an inverse Mellin 
transform of the moments. If enough moments are known (and can be continued to complex 
n), one can make a reliable reconstruction of the distributions numerically (see e.g. Ref. ||). 
However, since lattice data available at present (and in the foreseeable future) only exist 
for the first three nontrivial moments, this procedure becomes impractical. Alternatively, 
one can reconstruct the parameters in a given PDF parameterization by working with the 
moments directly. A similar method was applied to the reconstruction of polarized parton 
distributions in Ref. 0, while in Ref. |T(| the x dependence of PDFs was investigated subject 



to constraints on their small- and large-x behavior. In the present work, however, we wish 
to determine what constraints on the x dependence of the valence distribution are imposed 
by the actual lattice data. 



A typical parameterization of quark (or antiquark) distributions used in global fits [|TT^ [T3 
is of the form: 

xq(x) =ax b {l-x) c (l + ey/x + -fx) , (2) 
for which the corresponding moments are given by: 
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(x n ) q = a [B(l + c,b + n) + eB(l + c, § + b + n) + jB(l + c, 1 + b + n)] , (3) 

where B(x,y) is the /3-function. For a given set of PDF moments, the parameters a,b,c,e 
and 7 can be fitted using Eq. (|3]), and the corresponding x distribution obtained from Eq.(^). 
To test the effectiveness of this method, we calculate moments of PDFs from the recent MRS 
fTT| , CTEQ [jT2[ and GRV jOJ parameterizations, and examine how well the reconstructed 
distributions agree with the original parameterizations. 

In order to make best use of currently available lattice data, we concentrate on the 
nonsinglet valence distribution u v {x) — d v (x) = [u(x) — u(x)) — [d(x) — d(x)]. While the 
present analysis can be easily extended to the singlet sector, the singlet quark distribution 
on the lattice receives contributions from disconnected diagrams, corresponding to opera- 
tor insertions in quark lines which are disconnected (except through gluon lines) from the 
nucleon source. Disconnected diagrams are considerably more difficult to compute, so only 
exploratory studies of these have so far been completed fl4] . On the other hand, since the 



disconnected contributions are flavor independent (for equal u and d quark masses), they 
cancel exactly in the nonsinglet u — d combination. 

To enable comparison with the reconstructed distributions, we employ Eq. (0) to con- 
struct an 'average parameterization' from the global PDF analyses. Note that the PDF fits 
in Refs. UTTHT^I are given for u v (x) and d v (x) separately, whereas here we require a param- 



eterization of the difference. Using an equally weighted numerical fit to the u v (x) — d v (x) 
difference (at next-to-leading order in the MS scheme, at a scale Q 2 = 4 GeV 2 ) from the 
parameterizations of Ref. [|n|-|i~3|1, we find an average fit: 



x(u v (x) — d v (x)) aV g = a A x bA (l - x) CA (l + caVx + Iax) , (4) 

with 6a = 0.476, ca = 3.750, eA = —3.152, 7a = 16.64, and with oa = 0.561 determined 
by the normalization condition, Jq 1 dx (u v (x) — d v (x)) = 1. Figure [l] shows the average dis- 
tribution (long-dashed line), with the spread between the three parameterizations indicated 
by the shaded region. The lowest six moments of this average valence u v — d v distribution 
are calculated to be: 



(xXU = 1 , ^)uU = 0-163 , (xTZU = 0.0543 , 
(xTZ S - dv = 0-0229 , (xTZ g - dv = 0.0111 , (x 5 )Z S - dv = 0-00599 . 

We have investigated fits to these moments using the following seven fitting functions: (i) 
the full form of Eq. (Q) with unconstrained parameters; (ii) a simplified form with 7 = 0; 
(iii) a form with e = 0; (iv) both e = and 7 = 0; (v) with the constraint e = (vi) 
7 = 7a; and finally (vii) e = e& and 7 = 7a- Fits (iii) and (iv) are unable to reproduce the 
small x behavior of the distribution and are therefore discarded. The parameters determined 
by x 2 fits with each of the remaining forms give rise to similar distributions. We use the 
differences in the values of the parameters found from the various fits as an indication of the 
systematic error in the reconstruction procedure. In particular, for the parameter b (which 
is of special interest because of its relation to Regge phenomenology-see Section IV below), 
we find the best fit value b = 0.4 ± 0.2. 

Figure [l] shows the reconstructed x(u v (x) — d v (x)) distribution (short-dashed line) which 
is generated by the fit form (i) using all six moments above. With just the lowest 4 moments, 
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however (as are available from lattice data), this form does not produce an accurate fit, as 
one is attempting to constrain more parameters than the number of data points permitQ. 
This is illustrated by the dotted line in the Fig. [I]. In contrast, with fit (vii), which has only 
3 free parameters, one can obtain a fairly reliable reconstruction even when fitting to only 
the n = 0, 1, 2 moments, as indicated by the solid curve in Figure [I| 

As well as understanding the systematic uncertainty associated with the functional form 
of the distribution to be reconstructed, one must also quantify the error arising from the use 
of only a fixed number of moments. When the number of moments used with fitting form 
(vii) is increased to 4, 5 or 6 the fit improves by only a few percent. With four moments, 
fits (ii), (v) and (vi) also give acceptable results, similar to those shown for fit (vii). For 
simplicity, however, we will concentrate on fit form (vii) hereafter. 

The overall agreement between the original and reconstructed distributions is clearly 
excellent, and comparable to the difference between different parameterizations (shaded 
region in Fig. [I]). Having confidence that a reliable reconstruction is possible from just the 
lowest few moments, we next ask what do extrapolations of lattice data on the n = 0, 1, 2, 3 
moments tell us about the x dependence of the distribution u v (x) — d v (x)l 

III. QUARK DISTRIBUTIONS FROM LATTICE MOMENTS 

A crucial question for relating lattice calculations to phenomenology is the reliability of 
the extrapolation from the region of large quark masses (m q ) to the physical value. Typically, 
lattice data have been analyzed assuming a linear dependence on the quark mass. However, 
as pointed out in Ref. ||, if one is to extrapolate the lattice moments into a region where 
pion loops play an important role, it is vital to use an extrapolation formula which has the 
correct chiral behavior. In this section we first motivate an extrapolation formula which 
ensures both the correct chiral and heavy quark limits. Following this, we use the values of 
the moments at physical quark masses extracted from the lattice data to reconstruct the x 
dependence of the nonsinglet distribution. 

A. Constraints from the Chiral and Heavy Quark Limits 

The spontaneous breaking of the chiral SUk(Ny) xSUr(N/) symmetry of QCD generates 
the nearly massless Goldstone bosons (pions), whose importance in hadron structure is well 
documented. At small pion masses, hadronic observables can be systematically expanded 
in a series in m n |J|. The expansion coefficients are generally free parameters, determined 
from phenomenology or obtained from low energy models. One of the unique consequences of 
pion loops, however, is the appearance of so-called chiral logs. From the Gell-Mann-Oakes- 
Renner relation one finds that m% ~ m q at small m n , so that these terms are nonanalytic 
functions of the quark mass. Because the nonanalytic terms arise from the infra-red behavior 



1 Here we count the zeroth moment as a fit point with zero error, and include the overall normal- 
ization parameter a in our parameter counting. 
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of the chiral loops, they are generally model independent. For some observables these 
nonanalytic terms give rise to large deviations from linearity near the chiral limit [|l(| . 

The leading order (in m n ) nonanalytic term in the expansion of the moments of PDFs 
was shown in Ref. [|J to have the generic behavior m^logm,,-. In Ref. || a low order 
chiral expansion for the moments of the nonsinglet u — d distribution was developed which 
incorporated the leading nonanalytic (LNA) behavior of the moments as a function of m q . 
Here we extend this in order to connect also with the heavy quark limit, in which the valence 
quark distributions become (^-functions centered at x — 1/3: 

u(x) — d(x) m - — >° 5(x — |) . (5) 
The corresponding moments in the heavy quark limit are therefore given by: 

Wu-* m ^ ^ • (6) 

It is straightforward to incorporate this constraint into the extrapolation formula, and its 
inclusion makes no significant difference to the central results of Ref. |J, as illustrated in 
Fig. |2|. For completeness, we give the extrapolation formula which explicitly satisfies both 
the heavy quark and chiral limits: 

( jjl ^ \ jji ^ 

1 + c LNA ml log f- + b n n , (7) 



where the coefficient clna = _ (1 + 3g\) / (4:71 f n ) 2 Hl7| ) and b n is a constant constrained by 
Eq.(§: 

K = - a n (l - // 2 clna) ■ (8) 

In this extrapolation scheme there are therefore three free parameters for each moment: 
an, m b,n, and /i. We find that the dependence of the mass parameter m hn on n is quite 
small, and therefore fix mb, n = mb for all n. As the b n term in Eq. (|7|) is included in order 
to provide a linear dependence on m 2 in the region where lattice data exist, mb should be 
large enough such that this is indeed the case. For mb > 4 GeV the fits to the moments are 
essentially independent of m&; we therefore fix m& = 5 GeV. 

The mass fi essentially determines the scale above which pion loops no longer yield rapid 
variation. Analyses of a wide range of observables such as masses, magnetic moments and 
charge radii []TB[ as well as structure function moments suggest a common scale, /i ~ 0.5 GeV, 



which simply reflects the scale at which the Compton wavelength of the pion becomes 
comparable to the size of the hadron (without its pion cloud). We use the value [i = 550 MeV 
found in our previous fits to the lattice moments of the u — d distributions in Ref. ||. 
Consequently the extrapolation for each moment is essentially parameterized by only one 
variable. Nevertheless, the best fit using Eq.(^) is as good as that found in our earlier work 
H (which did not impose the additional constraint of the heavy quark limit), as illustrated 
in Fig. |2] — the latter is given by the short-dashed curve and the former by the solid curve. 
For details of the fitting procedure and the error analysis we refer to Ref. ||. 

Having motivated the functional form of the extrapolation formula to connect with both 
the chiral and heavy quark limits, we now use the extracted moments to calculate the x 
dependence of the nonsinglet valence u — d distribution from the lattice data. 
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B. Results 



Before extracting information on the x dependence of a particular PDF from the lattice 
moments, one must note that the crossing symmetry properties of the spin-averaged quark 
distributions (see Eq. ([!])) mean that the moments of u — d calculated on the lattice corre- 
spond to the valence u v — d v distribution only for even n. The odd n moments, on the other 
hand, correspond to the combination u + u — d — d = u v — d v — 2(d — u). The difference 
between the two represents the famous violation of the Gottfried sum rule, when d ^ u [|T],0 . 
Some care must be taken, therefore, if one attempts to use information on both even and odd 
moments to extract a valence quark distribution. To minimize possible contamination of the 
extracted valence distribution arising from the flavor asymmetry of the sea, in the present 
analysis we subtract the values of the phenomenological moments of the d — u difference, as 
given by global parameterizations of data [pT|-[T3|, from the calculated odd moments. 

Phenomenologically, the difference d(x) — u(x) is found to be concentrated at small x 
(d — u ~ (1 — a;) 14 at Q 2 ~ 50 GeV 2 from the Fermilab Drell-Yan data 0), and will therefore 
become progressively less important for larger n compared with the valence distribution. 
The largest effect is for the n = 1 moment and results in a correction of < 10% of the 
value of (x) u -d in Eq. (§), if one adds the phenomenological value 2 f Q x(d(x) — u{x))dx = 
0.017(2) P"T|-|T3| to the extrapolated first moment. For the n = 3 moment, the corresponding 
correction is a fraction of a percent. In principle, given sufficiently many moments, one can 
fit both the total u — d and valence u v — d v distributions separately from the n even and n odd 
moments, respectively — however, given the small number of moments currently available, 
the extraction of the valence distribution is the most feasible solution. 

Table I lists the resulting extrapolated values for the lowest three nontrivial moments of 
the total u — d distribution at the physical quark mass for both a linear extrapolation and 
the chirally symmetric extrapolation of Eq. (|7|). Also listed are the corresponding values of 
the valence fit parameters a, b and c from \ 2 fits. As discussed in Section II, because there 
are only four data points with which to fix parameters, we use Eq. (|3]) with the fitting form 
(vii), namely fixing e = e& and 7 = 7a- With these parameters, the resulting fits to the 
moments of the valence u v — d v distribution are displayed in Fig. [5] for both the linear and 
improved, chirally symmetric extrapolations. For clarity we plot 3™ times the moments, so 
that the horizontal line at unity represents the heavy quark limit of the moments. The lightly 
shaded region around the moments from the improved chiral extrapolation corresponds to 
a 1 standard deviation variation of the fit parameters from their optimal values. 

The corresponding distributions x(u v (x) — d v (x)) are displayed in Fig. |. Once again, 
the lightly shaded region represents a la deviation from the central values, while the darker 
band illustrates the spread between the global parameterizations from Refs. [P"T|-|i"3||. A 
comparison of the distribution reconstructed from the lattice data using the improved chiral 
extrapolation with the phenomenological distributions shows reasonably good agreement. 
On the other hand, the linear extrapolation gives a distribution (scaled by a factor 1/2 
in the figure) which is much smaller than the phenomenological distributions at small x, 
and consequently has a significantly higher peak in the intermediate x region, centered at 
x ~ 1/3 - reminiscent of a heavy, constituent quark-like distribution. 

The improved extrapolation formula (0) enables one to study not only the quark distri- 
bution at the physical quark mass, but also the dependence of the distribution on the quark 
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(or equivalently, pion) mass. In particular, it allows one to trace how the x dependence 
changes as one goes from the chiral limit to the heavy quark limit, where the distribution 
approaches a ^-function. Table [Tl] lists the extrapolated values of the total moments at 
various pion masses. Fitting the corresponding valence moments with the functional form 
(vii) from Section II, this table also shows how the resulting fit parameters change as the 
pion mass is increased from m n = to = 5 GeV. Notice that the fit parameters b and c 
become larger with increasing m n , as the x dependence of the valence distribution becomes 
less singular at small x. This is dramatically illustrated in Fig. |5], where the correspond- 
ing valence x{u v — d v ) distribution is shown in the chiral limit (m^ = 0), at the physical 
mass (m„- = 0.139 GeV), and at = 0.5, 1 and 5 GeV. As increases, the distribution 
becomes visibly more sharply peaked, with the peak moving towards the limiting value of 
x = 1/3. 

Quite remarkably, the distribution at = 5 GeV resembles the distribution in Fig. |] 
(dot-dashed curve) extracted from the linearly extrapolated moments. Since the linearly 
extrapolated moments are not constrained by the expected behavior in the heavy quark 
limit, it is noteworthy that a constituent quark-like distribution, peaking near x ~ 1/3, 
nevertheless arises if the moments are flat as a function of m n , as one would expect within 
a constituent quark model. This in fact provides further a posteriori justification for the 
choice made in Eq.(|7]) to build in the correct heavy quark limit (as well as the correct chiral 
behavior). 



IV. REGGE BEHAVIOR AND THE J PC = 1 TRAJECTORY 



In the previous section we extracted the exponent b which governs the small-x behavior 
of the nonsinglet valence quark distribution u v — d v . According to Regge theory [[Tj|, this 



exponent is related to the slope of the isovector J PC = 1 trajectory, corresponding to the 
exchange of p (1 ), a<i (2 ++ ), p 3 (3 ) and so on. If one makes the reasonable hypothesis 
that Regge theory also holds for quarks with masses in the range m q ~ 30-190 MeV used 
in lattice simulations, one can estimate the masses of mesons on this trajectory from fits to 
the valence quark moments versus m q . 

From Regge theory the elastic scattering amplitude at high center of mass energy, s, 
behaves as: 

A(s,t)~s a ^, (9) 

where a(t) is a function describing the meson trajectory: 

a(t) = a + a' t , (10) 

with a' the slope and «o the intercept of the trajectory. The optical theorem then relates 
the imaginary part of the scattering amplitude at forward angles, A(s,t = 0), to the total 
cross section, a: 

a ~ ^ ~ s^ 1 . (11) 

At large s the approximate energy independence of total cross sections is usually attributed 
to the exchange of the Pomeron (IP) trajectory, which has intercept a^" ~ 1. On the other 
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hand, the isovector cross sections (corresponding to nonsinglet distributions) are dominated 
by the exchange of a trajectory with qq quantum numbers corresponding to the p meson 
and its excitations, which has intercept «q ~ 0.5. 

In deep-inelastic scattering at large s 3> Q 2 , corresponding to x — Q 2 /2p-q ~ Q 2 /s <C 1, 
the total virtual photon-nucleon cross section is dominated by the singlet quark and gluon 
distributions, and behaves as a 7 *^ ~ ~ x 1_Q o . The nonsinglet distribution, q^s, on the 
other hand, is governed by the exchange of the J PC = 1 trajectory, xq^s ~ x 1 ~ a o ~ x 1 ! 2 
(recall that the best fit to the average u v — d v difference in Eq. (|4]) gave an exponent b = 
0.476). Identifying ctg ~ 1 — b from the previous section as a function of quark mass, one 
can therefore extract information on the m q dependence of the Regge intercept from the m q 
dependence of the nonsinglet quark distribution. 

Having determined the intercept b « 1 — «q from the PDF fits, one needs to determine 
the slope of the trajectory as a function of m q . In the infinite mass limit, orbital excitations 
of mesons become energetically degenerate with the L = state. Within the Regge picture, 
this is possible only if the Regge intercept «q - — ► — oo, or b - — > oo, which is consistent 
with the valence distribution approaching a 5-function. One expects, therefore, that the 
slope a' should increase as m q increases from its physical value. The behavior a(J - — > — oo 
is indeed what is seen from the dependence of b in Table II. The slope at larger, but 
finite, values of m q is most accurately determined from the mass of the p meson, which has 
considerably smaller errors than the masses of higher lying mesons on the same trajectory. 
Using the lattice data for the p meson and the values of a.Q{m q ) generated from the fits 
in Section III, one can then make predictions for the behavior of the masses of the orbital 
excitations function of m„. 

In Fig. H we show the predicted 1 Regge trajectory at m n = 0.785 GeV (solid line), 
compared with the trajectory at the physical light quark mass (dashed). A fit through the 
central values of the parameter b and the p mass at m n = 0.785 GeV (as calculated in 
Ref. |TP|]) yields a slope which is somewhat larger than that of the trajectory at the physical 



quark mass, which is consistent with the expected trend towards the heavy quark limit. 
The darker shaded region represents the statistical error in the extrapolation of the lattice 
moments, obtained from fitting the upper and lower bounds of the original lattice data. 
The lighter region indicates an estimate of the systematic error associated with the fitting 
procedure as discussed in Section II. 

Although lattice data for the masses of orbital excitations are scarce, there have been 
some pioneering calculations of masses of the a 2 (J PC = 2 ++ ) and p 3 (J PC = 3 ) mesons 
by the UKQCD Collaboration ||19|| , indicated by the filled boxes in Fig. || Comparing with 



the predictions from the PDF analysis, the lattice data are in reasonable agreement: the 
calculated p3 meson mass lies within the predicted band, albeit within large errors, while 
the ci2 mass lies on the edge of the predicted range. Needless to say, further exploration 
of the masses of excited mesons within lattice QCD would be very helpful in testing these 
predictions. 

V. CONCLUSION 

In this paper we have extracted the x dependence of the unpolarized valence quark 
distribution x(u v (x)—d v (x)) from lattice QCD data on the moments of the u — d distribution. 
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After establishing that the basic features of the x distributions can be reconstructed from 
just the lowest four moments, we analyze the lattice data using an extrapolation formula 
which embodies the correct behavior in both the chiral and heavy quark limits. 

The importance of ensuring the correct chiral behavior is illustrated by comparing the x 
distributions obtained by extrapolating the lattice data using a linear and an improved chiral 
extrapolation. While the improved extrapolation gives an x distribution which is in quite 
good agreement with the phenomenological parameterizations, the linearly extrapolated data 
give rise to distributions with the wrong small-a; behavior, which translates into a much more 
pronounced peak at x ~ 1/3, reminiscent of a heavy, constituent quark-like distribution. 
The approach to the heavy quark limit is explicitly mapped out by studying the dependence 
of the x distributions on the pion mass, from the chiral limit to pion masses of several GeV. 

Our analysis suggests an intriguing connection between the small-a; behavior of the va- 
lence distributions and the m q dependence of meson masses on the J PC = 1 trajectory. 
While the existing lattice data are in qualitative agreement with the predictions from the 
PDF analysis, future results on the masses of excited meson states at varying quark mass 
will enable the accuracy of this relation to be tested more thoroughly. Plans by a number 
of lattice groups to study the spectrum of excited hadrons should provide valuable insights 
into this problem. 
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TABLES 



Extrapolation 


(x 1 ) 


<x 2 ) 


<x 3 > 


a 


b 


c 


Linear 


0.262 


0.0843 


0.0340 


658 


3.51 


10.5 


Improved 


0.149 


0.0461 


0.0200 


0.439 


0.351 


3.77 



TABLE I. Moments of the total u — d distribution obtained from a linear extrapolation and 
the improved chirally symmetric extrapolation, Eq. (^), to the physical quark mass. Also shown 
are the best fit parameters a, b and c for the valence u v (x) — d v (x) distribution, using Eq. (||) with 
e = eA and 7 = 7a constrained to the average values (see text). 



m w (GeV) 


(x 1 ) 


<x 2 > 


(x 3 ) 


a 


b 


c 





0.121 


0.0374 


0.0161 


0.282 


0.245 


3.55 


m phys 


0.149 


0.0461 


0.0200 


0.439 


0.351 


3.77 


0.5 


0.223 


0.0689 


0.0296 


2.39 


0.957 


4.96 


0.785 


0.248 


0.0769 


0.0328 


7.30 


1.45 


5.91 


1 


0.258 


0.0801 


0.0341 


13.8 


1.74 


6.48 


5 


0.303 


0.0978 


0.0367 


125000 


6.20 


15.4 



TABLE II. Moments of the total u — d distribution obtained from the improved chiral extrap- 
olation formula, Eq. ([?]), at various pion masses; and best fit parameters for the corresponding 
valence u v {x) — d v {x) distribution using Eq. (||) with e = €a and 7 = 7a constrained to the average 
values. 
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FIGURES 




x 



FIG. 1. Quality of reconstruction of the valence x(u v {x) — d v (x)) distribution from several low 
moments: the shaded region represents the spread between different next-to-leading order distri- 
butions from global parameterizations |Tl|-^| (at Q 2 = 4 GeV 2 in the MS scheme), while the 
long-dashed line represents a parameterization of the average of the three distributions, Eq. (H). 
The short-dashed line (which is almost indistinguishable from the long-dashed, average parame- 
terization) is the distribution reconstructed from the lowest six moments of the average parameter- 
ization using Eq. (||) with e and 7 unconstrained. The dotted curve indicates the fit obtained when 
only four moments are used with the same fitting form. In contrast, the solid lines represent the 
distribution reconstructed from the lowest three moments (n = 0, 1, 2) using Eq. (||) with e and 7 
constrained to the values obtained from direct fits to the average distribution, e = eA and 7 = 7a- 
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Experiment 




FIG. 2. Moments of the u — d distribution. The lattice data [|6|,|7j] are fitted using a linear 
function (long-dashed), and with the improved chiral extrapolation in Eq. ([?]) (solid). Also shown 
is the fit from Ref. || (short-dashed) which does not impose the heavy quark limit, Eq.(^). 
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FIG. 3. The lowest four moments of the valence u v — d v distribution (scaled by 3 n ) at the 
physical quark mass, extracted from the fit to the lattice data using a linear extrapolation (di- 
amonds) and the improved chiral extrapolation, Eq. (Q) (stars). The solid and dot-dashed lines 
are x 2 fits to the improved and linearly extrapolated moments respectively using Eq. (||), with e 
and 7 constrained to their average values. The shaded region represents a 1 standard deviation of 
the fit parameters about the optimal values for the improved extrapolation. The short-dashed line 
represents the heavy quark limits of the moments. 
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FIG. 4. The physical valence x(u v (x) — d v (x)) distribution, extracted using the improved chiral 
extrapolation of the lattice moments (solid), and a linear extrapolation, scaled by a factor 1/2 
(dot-dashed). The lighter shaded region indicates a la variation of the fit parameters about the 
optimal values for the improved extrapolation, while the dark shaded region represents the spread 
between global parameterizations [lT|-I^|. 
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FIG. 5. The nonsinglet valence distribution x(u v (x) — d v (x)) extracted from the improved 
extrapolation formula, Eq. (0), for various pion masses: = (short-dashed), m n = 0.139 GeV 
(solid), m w = 0.5, 1 and 5 GeV (long-dashed). The fit parameters are tabulated in Table II. 
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FIG. 6. Regge plot of the spin, J = a(j, versus t = (mass) 2 of mesons on the J PC = 1 
trajectory, at the physical pion mass (dashed) and at m n = 785 MeV (solid), with the point at 
t = obtained from the best fit to b (see Table II). The darker shaded region, which includes the 
error on b, illustrates the statistical error in the extrapolation of the lattice moments, while the 
lighter shader region indicates the systematic error in the fitting procedure. The physical masses 
of the p, a% and p$ mesons are indicated by stars, while the boxes represent lattice masses at 
m n = 0.785 GeV from the UKQCD Collaboration @. 
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